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There is increasing evidence that conventional cold dark matter (CDM) models lead to conflicts 

between observations and numerical simulations of dark matter halos on sub-galactic scales. Spergel 

and Steinhardt showed that if the CDM is strongly self-interacting, then the conflicts disappear. 

However, the assumption of strong self-interaction would rule out the favored candidates for CDM, 

namely weakly interacting massive particles (WIMPs), such as the neutralino. In this paper we 

propose a mechanism of non-thermal production of WIMPs and study its implications on the power 

("^ ' spectrum. We find that the non-vanishing velocity of the WIMPs suppresses the power spectrum 

f^ I on small scales compared to what it obtained in the conventional CDM model. Our results show 

^^ ■ that, in this context, WIMPs as candidates for dark matter can work well both on large scales and 

Cn ' on sub-galactic scales. 

Q There is strong evidence for the existence of a substantial amount of cold dark matter (CDM) . The leading candidates 

D for CDM are weakly interacting massive particles (WIMPs), such as the neutralino. The neutralino is the lightest 

II ] supersymmetric particle. In models with R parity it is stable, and its mass density in the universe is generally assumed 

. ■ to be a relic of an initially thermal distribution in the hot early universe. Assuminig, in addition, the presence of a 

small cosmological constant, the CDM scenario is consistent with both the observations of the large scale structure 

fT^ . of the universe (^ IMpc) and the fluctuations of the cosmic microwave background [Q. Many experiments searching 

J> ' for neutralino dark matter particles are under way. 

CO , The coUisionless CDM scenario, however, predicts too much power on small scales, such as a large excess of dwarf 

^^ ■ galaxies [^,|| , the over-concentration of dark matter in dwarf galaxies [§-g and in large galaxies . Recently Spergel 
and Steinhardt proposed a new concept of dark matter with strong self interaction M . This puts WIMPs as candidates 
for dark matter in considerable jeopardy P,[lO|. 
p^ ' In this paper we propose a scenario with non-thermal production of WIMPs. These WIMPs could be relativistic 

f^ ■ when generated. Their comoving free-streaming scales could be as large as of the order 0.1 Mpc or larger. The 
density fluctuations on scales less than the free-streaming scale would then be severely suppressed. Consequently 
the discrepancies between the observations of dark matter halos on the sub-galactic scales and the predictions of the 
r^' standard WIMPs dark matter picture could be resolved. 

P I To begin with, we consider a general case of non-thermal production of the neutralinos by the decay of topological 

defects such as cosmic string fl^, by the decay of an unstable heavy particle, or produced non-thermally by the 
^ [ reheating process in a scenario of inflation at low energy scale ||l^ (see also [|3|). The momentum distribution 
• • function of the neutralinos is for simplicity assumed to be Gaussian: 

where Pc is the central value and a describes the width of the distribution. 

Given a model, the parameters Pc and a can be determined. For instance, in the supersymmetric version of the 
Ub-l{^) model, the (Higgsino-like) neutralinos arise directly from the decay of the right-handed neutrinos and their 
superpartners uM. In this model, Pc is about a half of the mass of the mother particles. For a two-body decay, if 
the mother particle is at rest, the distribution function f(p) is a J- function. The value of a characterizes the average 
non- vanishing velocity of the mother particles (when a -^ 0, f{p) approaches a (5- function). In the model of Ref. [O, 
since the mother particles when released from the cosmic string loop are non-relativistic, a is small compared with pc 
and one has Pc ~< p >~< p^ >^/^. 

In Eq. (|0.l| ) , A is a normalization factor determined by the energy density of the non-thermal component 

PNT^^^ j E{p)f{p)p'dp, (0.2) 

where E{p) = (p^+m^y^^ and m is the rest mass of the dark matter particle. Given that the physical momentum p{t) 
scales as the inverse of the cosmic scale factor a(i), we define r = a{t)p{t)/m. During cosmic evolution r is a constant. 
Throughout this paper we set a(io) = 1, so r can be understood as the velocity of the particles at the present time 
(note that the dark matter particles are non-relativistic now even though they are relativistic when generated). 
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The comoving free-streaming scale Rf for the non-thermal particles can be calculated as follows ||TJ,|l^: 



Rf = / ^hrdt' ^ / ^hrdf 



" ^--^^^(1 + ^^^r In [^1 + ^g(^ ;^^^), + ^^(T^J > (0.3) 

where Tc = a(t)pc(t)/m — Pc{to)/m and the subscript '^q' denotes radiation-matter equality. Below the free-streaming 
scale, the power spectrum will be severely damped. To account for the lack of substructure in the Local Group, N-body 
simulations study show that the free-streaming scale of the dark matter should be ~ O.lMpc |l^. This corresponds 
to Vc ~ 10~^, which gives rise to a constraint on the parameters of our model. 

To calculate the power spectra for our models, we adapt the CMBFAST code [Q so that it applies to a thermal as 
well as a non-thermal dark matter distribution. To check the program we replaced the non-thermal distribution func- 
tion in our adapted code by the standard thermal distributions for hot dark matter and cold dark matter respectively, 
and calculated the power spectra, obtaining the same results as with the standard CMBFAST code. 

In the numerical calculations, we took the presently favored values ilm = 0.4, J7a — 0.6, h = 0.65 and flbh^ — 0.02 
of the relevant cosmological parameters, where flm, ^A; and f2f, are the ratios of the contributions from total mass, 
vacuum energy, and baryons to the total density of the universe. In general, there are two sources of WIMPs which 
contribute to firm one is from thermal production, the other is from non-thermal production. For simplicity, in our 
study we neglected the contribution of WIMPs from the thermal production mechanism. 

In Fig.l, we show the power spectrum of our model with Vc = 1.5 x 10~^. For comparison, we also plot the power 
spectra for the conventional CDM model and the warm dark matter (WDM) model with mw — IkeV (For this 
mass the free-streaming scale is O.llMpc ||l^). Throughout this paper, we normalize all power spectra by the COBE 
observations [[19[ . 

We have varied the parameters of the non-thermally produced WIMPs dark matter model in the numerical calcu- 
lations. Different values of re gives rise to different power spectra. However the spectrum is independent of the rest 
mass m of the dark matter particle (for fixed Tc), and is very insensitive to a for a <$^ pc. This can be understand 
easily from Eq.(3.3) which shows that the comoving free-streaming scale Rf only depends on r^. 



One can see from Fig. 1 that on large scales, the three power spectra are the same, which shows that our model 
retains all the merits of the conventional CDM model on large scales, while on sub-galactic scales, the power spectrum 
of non-thermal dark matter (NTDM) model is damped severely relative to that of the CDM model, and is close to 
that of the WDM model with mw — IkeV. Ref. |1^ argued that such a warm dark matter scenario provides a solution 
to the substructure problem. 

The model we present in this paper with non-thermal production of WIMPs provides a promising scenario for 
large-scale structure formation of the universe. However, we need to check its consistency with observations on small 
scales, especially on scales of the Lyman-a forest. In Fig. 2 we plot the power spectrum of our model and the observed 
power spectrum of the Lyman-a system at z = 2.5 shown as the filled circles with error bars pO|j2l[ |. For comparison 
we also give the power spectra for the conventional CDM and WDM models. In fitting to the observed data, we 
choose the primordial spectral index n — 0.97 for all models. The mass for the WDM particles is choosen as 750eV, 
and the parameters for the NTDM models are re = (1.3, 1.4, 1.5) x 10~^, respectively. 

To obtain a quantitative constraint on Vc, we closely follow the study by Narayanan et al in Ref. ||l^. They studied 
explicitly the power spectrum of WDM models and came to the conclusion that any mechanism which suppresses 
the conventional CDM linear power spectrum more severely than a 750eV WDM particle will be inconsistent with 
Lyman-a forest observations in the scale range —2.5 < lg[k(km~^ s)] < —1. This requires that the linear theory power 
spectrum of any substitute to the conventional CDM models be higher than that of a WDM model with mw — 750eV 
p8[ . For a flat universe with f^A = 0.6 the range of these scales corresponds to 0.4ft,~^ Mpc< k < 12.8 h^^ Mpc. 
From Fig. 2 one can see that the larger the value of Tc is, the lower the small-scale power spectrum becomes. By 
comparing the values of the power spectra at the upper limit of the above range in k with the power spectrum of the 
WDM model with mw = 750eV, we obtain an upper limit on Vc of Vc < 1.5 x 10~^. 

Next we examine the constraints on our model by the recent studies of the phase-space density |2^-Q. The 
particle phase-space density Q is defined as Q = p/ < v^ >^/^, where p is the energy density and < u^ > is the mean 
square value of the particle velocity. The astronomically observable quantity is the mean coarse-grained phase-space 
density. In the absence of dissipation, the coarse-grained phase space density can only decrease from its primordial 
value. Thus, one can use the observed maximum phase-space density to set a lower limit on the phase-space density 
for the dark matter particles. The highest observed phase-space density is obtained from dwarf spheroidal galaxies: 
Qobs ^ lO~''M0pc^(km/s)~^ |Q. For our models, < v"^ >~ r^ at the present time, and therefore 

Qo ~ PNTo/rl . (0.4) 



Because the primordial phase space density decreases with time when the particles are relativistic and becomes a 
constant after the particles become non- relativistic, one requires Qo > Qobs, which can be translated into a constraint 
on Vc'. 

re < {pNTo/Qobsf^ - 2.5 X 10-" . (0.5) 

This limit is of the same order of magnitude but slightly weaker than that from the observations of the Lyman-a 
forest. 

In conclusion, we have demonstrated that our model with non-thermally produced WIMPs is consistent with 
observations on scales ranging from the Lyman-a systems to the CMB. Compared to the conventional CDM model, 
the power spectrum of our model on sub-galatic scales is severely damped, which will account both for the lack 
of substructure in the Local Group and for the observed smooth inner halos[J. In addition, during the structure 
formation, the non- vanishing velocity of the WIMPs may further reconcile the discrepancies between the theoretical 
predictions and observations on sub-galactic scales |^,Q. However, the detailed effect of the dark matter velocity 
on the formation of halos needs further study with high resolution N-body simulations. 

To solve the sub-galactic-scale problem, our study in this paper requires the value Vc ^ 1.5 x 10^. This has 
implications on the models of non-thermal production of WIMPs. In the scenario of Ref. [O where neutralinos 
are released by cosmic string decay, Vc can be approximately expressed as Vc — {TQ/2m){M/Td), where Tq is the 
present temperature, Td is the temperature when the string decay first leads to nonthermal neutralinos, and M is the 
characteristic mass of the mother particle and the scale of the defect which determine the initial momentum. For m 
in the range of 50Gev to lOOGev, M is required to be in the range of 1.6 x lO^Gev to 6.4 x lO^Gev. This is consistent 
with the results given in Ref. |ll[ and shows that the value which we used in the text is quite realistic. 

In summary, the discrepancies between theory and observations on sub-galactic scales disfavores the conventional 
WIMPs CDM model. In this paper we show that if the dark matter particles have a non-thermal origin, these 
discrepancies can be resolved and the WIMPs remain good candidates for the dark matter particle. 
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FIG. 1. Comparison of the power spectra of the CDM model(long dashed curve), the WDM model with mw ~ lkeV(short 
dashed curve) and the NTDM model with re = 1.5 x 10^^ (solid curve). 
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FIG. 2. The power spectra of the CDM model (long dashed curve), the WDM model with mw = 750eV (short dashed curve) 
and the NTDM models with Tc — (1.3, 1.4, 1.5) x 10^^ (solid curves, from top down), compared to the observed lyman-a P{k) 
&t z — 2.5 (filled circles with error bars). 



